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The common junction of three dissimilar general shells of revolution is analyzed. Axisym- 
raetric loading may be in the form of surface forces, concentrated forces and moments at the 
junction, and arbitrary thermal gradients. Basic differential equations available for elastic 
shells are extended for application in the elastic-plastic regime. The von Mises-Hencky 
yield criterion, deformation theory of plasticity, and successive approximations are used to 
determine plastic strains. Postyield material behavior is arbitrary. Linearized finite-differ- 
ence equations are solved directly using coefficient matrices that incorporate conditions of 
equilibrium and compatibility at the junction discontinuity. The solution to a sample prob- 
lem is given. 


Nomenclature 

= extensional rigidity, f h E d£, lb/in. 


= flexural rigidity, 


(1 


■h >j> w ' 


lb-in. 


«o 

0 

A( 


= modulus of elasticity, psi 
V = radial and tangential stress resultants, lh/in. 

= thickness of shell wall, in. 

= couple, lb-in. /in. 

— distributed surface load, psi 

z = radial, meridional, and axial coordinates to reference 
surface, in. 

= temperature, °F 

— radial deflection, in. 

*= axial stress resultant, — — | rcmpvd£, lb/in. 

= coefficient of thermal expansion, in. /in .- °F 
= reference length, in. 

— meridional rotation, rad 
) = increment of ( ) 

— normal strain, in. /in. 

=* transverse coordinate, satisfies J* h dt; = 0, in. 

— change in curvature, (in.)” 1 
= Poisson’s ratio 

= «/«• 

= normal stress, psi 


Subscripts 

a,b,c = junction, shells a, 6, and c, respectively 

d = boundary, shell c 

e,t t p = effective, total, and plastic 

H,V ~ radial and axial 

i y j = station indices 

— circumferential, meridional, and transverse 

Superscripts 

e = external 

(A) = reference surface 

Introduction 

I N certain missile and space vehicle structures, small plastic 
strains may be permitted in an effort to use material more 
efficiently. In addition to increasing allowable loads on 
continuous shell regions, a small amount of yielding may have 
the beneficial effect of relieving maximum stresses in regions 
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of discontinuity. The purpose of this paper is to present a 
shell analysis with which both this reduction in stress con- 
centration and the corresponding increase in strain concen- 
tration may be determined quickly and accurately. 

Many analyses have been presented for elastic shells of 
revolution containing discontinuities. One of the most gen- 
eral and convenient is that proposed by Radkowski et a!. 1 
In this reference, basic differential equations derived by 
Reissner 2 are solved by a very direct finite-difference tech- 
nique, which employs coefficient matrices. Continuous shell 
regions are analyzed separately for five independent sets of 
boundary conditions and are then joined at discontinuities 
using the well-known method of simultaneous junction equa- 
tions and superposition. 

Mendelson and Manson 3 presented a method for analyzing 
thermal stresses in a uniform cylindrical shell in the elastic- 
plastic regime. This method was later applied to general 
shells of revolution by Stern, 4 using the earlier Reissner- 
Meissner differential equations. Both of these references 
are applicable only to continuous shell regions. In addition, 
the numerical method employed in Ref. 4 requires the simul- 
taneous solution of a set of equations. Because each finite 
difference station along the shell meridian adds two equations 
to this set, this method is more limited than that of Ref. 1. 

References 1 and 2 were used by Wilson and Spier 5 to 
develop a solution for the small finite elastic deformations of 
continuous shells of revolution. This problem is similar 
to that of elastic-plastic deformation in that nonlinear terms 
occur in the basic differential equations of each. Also, in 
both Refs. 3 and 5 the equations are linearized by treating 
the nonlinear terms as knowns and correcting their values by 
successive approximations. 

The present analysis is much more general than previous 
work because it is applicable to both continuous and dis- 
continuous shell regions in the elastic-plastic regime. A 
wide range of shell geometries, loadings, and material prop- 
erties may be specified. For elastic discontinuity problems, 
multiple solutions and simultaneous junction equations are 
eliminated. In addition, a great measure of the generality, 
directness, and convenience of the shell analysis of Ref. 1 
is preserved in the presence of discontinuities and plastic flow. 

Method 

The common junction of three dissimilar general shells of 
revolution is used as a general discontinuity model (Fig. 1). 
Wall thicknesses, material properties, distributed loads, and 
temperatures may vary along the several shell meridians. 
In addition, material properties and temperatures may 
vary arbitrarily through the shell walls. These data may 
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Fig. 1 Meridional geometry of general shell discontinuity 
structure. 


be specified at discrete points, if desired. An external axi- 
symmetric force and couple are applied at the junction for 
further generality. 

The analysis of the three continuous shell regions is con- 
sidered first. Basic differential equations of Ref. 2 are ex- 
tended for application in the elastic-plastic regime in accord- 
ance with Ref. 3. These modified linearized equations are 
then solved by the finite-difference technique of Ref. 1 for 
two quantities, ( rH) and f3, from which all stresses and strains 
in the shell may be computed. The shell meridians are 
travereed twice: first, to compute coefficient matrices at 
each station using geometrical and load data; and, second, 
to calculate ( rH) and /?. 

After the continuous regions are analyzed, junction equa- 
tions are derived which express the structural continuity 
that exists at the intersection of the three shells. Since dis- 
continuity forces and couples are not linear with applied 
loads in the elastic-plastic regime, standard methods of 
analysis which use linear edge influence coefficients 1 * 7 are 
not applicable. Instead, the junction equations are satisfied 
by incorporating them directly into the previously mentioned 
coefficient matrices, which preserves the direct character 
of the finite-difference technique. Finally, the numerical 
technique is summarized, and a sample problem is solved for 
both strain-hardening and perfectly-plastic material behavior. 


v 




Fig. 2 Stress resultants, couples, and loads. 


In addition to the usual assumptions of thin shell theory, 
the von Mises-Hencky yield criterion and the deformation 
theory of plasticity are assumed to be valid. Postyield mate- 
rial behavior is arbitrary. 


Analysis 

A. Shell Equations 

2. Basic differential equations 

In accordance with the procedure developed by Reissner, 2 
equilibrium, compatibility, stress-strain, and strain-displace- 
ment equations in a thin elastic shell of revolution may be 
combined and reduced to two second-order coupled differen- 
tial equations. The dependent variables are a representative 
deformation ft and a representative stress function (rH). 
The independent variable is {, the meridional coordinate. 
If extended to include plastic as well as elastic and thermal ' 
deformations, these two basic equations become nonlinear 
and have the following form: 

(rHY + T(rH) f + 9 (rH) + A0 = \i + n + n (la) 

P" + W + *0 + MrH) = X 2 + r 2 + tt 2 (lb) 


Primes denote differentiation with respect to £, and the 
coefficients are 


r = 


(r/CaoY 

r/Ca 0 


T = 


(rD/apY 

rD/ao 


* # . z' 

V = 


A = - 


r/Ca 0 


Xl= _ir vpy +v t 

(L r/Cao r , 


rD/otQ 

(ratpu) + (raop H )'|+ 




x 2 — 


1 


[r'(rV)) 


rD/ao 

v) 


" - o^d{L e ^«) 

■ sk {[HW* Ei,y + df ] - 

i L v t f h E ^ p + | 


7Tl 


7T2 


(lc) 


The integrals through the shell walls may be evaluated 
analytically or numerically. 

The left-hand sides of Eqs. (la) and (lb) and the loading 
terms Xi and X 2 are identical with those derived in Ref. 2 for 
elastic shells. The thermal terms ti and r 2 are developed in 
Ref. 1 . The plasticity terms tti and 7 r 2 are derived in a similar 
manner, by using the elastic-plastic stress-strain relations of 
Ref. 3: 


<?$ = [E/( 1 - v 2 )][(«, + ve*) “ 

(eev + ve £p ) - (1 + v)aT\ 

<r £ = [E/( 1 - v 2 )][(e* + vee) - 

(<*p + ve ep ) - (1 + v)aT] 


(Id) 
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Although these last two terms contain nonlinear plastic 
strains, it is convenient to treat them not as unknown but as 
known quantities that may be determined in an iterative 
manner from previous, less accurate solutions. For example, 
Ti and 7 t 2 may be assumed to be equal to zero to obtain an 
elastic solution to Eqs. (1). Values of e$p and te v may be 
calculated from these elastic results and used to obtain some-" 
what more accurate approximations for ti and ir 2 . This 
procedure is repeated until additional solutions produce 
negligible changes in the plastic strain terms. 

General boundary conditions on Eqs. (1) can be expressed 
in terms of the quantities ( rH ), (rH) 0, and /?'. Because 
these equations will be solved using the method of finite 
differences, it is more convenient to express these boundary 
conditions directly in difference notation. This is done in the 
following section. 

2, Basic difference equations 

The shell reference meridians are divided into equal 
increments (Fig. 1 ) . Using central differences, the derivatives 
in Eqs. (1) may be expressed as 

( m , _ (rH)i+i ~ 2 (rfl),- + (r/Qj-i 
\rn)j - ,erc. 


Substituting Eqs. (4b) into (2a) gives 

( AjP ,• ~ 1 + Bj — F jPj_i)yj -f 2, a — 1 

(AjP j~ l qj + F — gi) = 0 j = a + 2, b — 1 

k c + 1, d — 1 

Tliis equation is satisfied independently of y } by setting 

-AiPr* + Bi - F jPj^i = Afr'qj + 

^ ' F &-i - 9i = 0 (4c) 

Therefore, 

Pi = (Bi - FjPj-i)~ l Aj (5a) 

f 2, o - 1 
j= a + 2,5-1 
U + i,d-i 

Qi = PjAj~ l {Qi — Ffti-. i) (5b) 

From Eqs. (3) and (4a), 

Pi = J qi = k P o+i = L q m+l = m (5c) 

Using the last of Eqs. (3) and Eq. (4a) results in 

y* = (/ - RPd-i)^(t - Rq^O (6) 


In matrix form, the basic shell difference equations then be- 
come 


( 2,o-l 

AjVi+i + Bj-yj + F iVj-i = gj j = < o + 2, 6 — 1 (2a) 

U + l, d- 1 


Capital letters represent 2X2 and lower-case letters 2X1 
matrices. 

These are 


A 

B 

F 

9 

y 


ri + (A£/2>r:' • o i p„ 01 
L 0 U + (A£/2)tJ Lo On J 

r-2 + (A*)*e (A©*A 1 pu 6u"| 

L (AflV y -2 + (A©*$J L&a 6»J 

fl ~ (A{/2)l>>- v 0 1 pii 0] 

L 0 1 —A£/2 ij LO fa] 

r(Ae*(Xi + n pq 

L(A^)*(Xj + Ti -f TT 2 )J L^J 



(2b) 


Boundary conditions that may involve first derivatives of 
(rH) and @ require the use of forward and backward differ- 
ences. A boundary condition may then be expressed as a 
general linear combination of the unknowns at the boundary 
and at the adjacent station. 

In matrix notation, the boundary conditions become 

?/i = k — Jy t 1 

2 /«+i = m — Ly m+i > (3) 

2/ d — t — Rpd-i j 


in which I equals the 2X2 identity matrix. 

Equations for P CJ q Ci yb , and y a are needed to complete this 
numerical integration. They will be derived in the next 
section from the junction conditions. With these additional 
equations, consecutive calculations of the y matrices may 
be made using Eqs. (6) and (4a), after the P and q matrices 
are calculated by Eqs. (5). For convenience, matrix equa- 
tions (3, 4a, 5, and 6) are expanded in Ref. 6. 


B. Junction Equations 
2. Derivation 

Certain conditions of compatibility and equilibrium must 
be f ulfill ed at the junction of the three shells to maintain 
structural continuity. These are 1) the junction is hinge- 
less; 2) the reference surfaces of the three shells experience 
equal radial deflections at their junction; and 3) a differen- 
tial ring element that includes the junction section is in static 
equilibrium. 

The first condition leads directly to the equations 

/*« = ft = ft (7a) 

If the superscript ( A ) is used to denote a reference surface, 
the radial deflection of this surface at any station is T&, j. 
Therefore, the second continuity condition may be expressed 
by the equations 


ed,a — e6,b — 60, e (7b) 

Referring to Fig. 2, the junction equilibrium conditions are 
(rH). + (rH)b - (rH). - fW (7c) 

in which f is the nominal radial coordinate to the junction, and 


3. Numerical integration 

Following the method of Ref. 1, the equation matrix (2a) 
is solved as follows. Set 


Vi = Qi - PjVi+1 

(4a) 

y,-i = g,_i - Pj-iyj ) 

(4b) 

y,+i = Pr'(9i ~ Vi) ) 



Mz,a + Mi, 6 = Mi, c - Mi (7d) 

Conditions (7b) and (7d) must be expressed in terms of 
(rH) and ft. Following closely the derivation in Ref. 1. 
modified slightly by the stress-strain equations (Id), 

« = A [— -- £ (rfl)"| + X, + r, + n 
C L °0 «0 r J 

t£ = -L r_„^ + lL (rfl)l + A, + 7-4 + *4 

* C L a. Oof J 


(8a) 
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junction equations become 



[ 

[ 


rp a - — (rF)l 

OoT J 

-vrpa + — (rV)l 
ot<>r J 


f t E„Td! 


U = TJ 


(r/0. = T~ [ SUrH ). + - S 2 ,„)/3 C + - 5,,.)] (10a) 

di.o 

(ri/)» = [5i,.(rtf). + («,.. - 5 2 , 6 )/3. + («,.. - «,.»)] (10b) 

di.b 

(rff). + (rH) b = (1 - f% x d l , c )(rH) e - 

f(fs u 5 2 ,c + s w )j3 c — f(fSu5 3tC + i/) (10c) 



k Etep ^ ^4 - -p- J h Ee$ p df 

and 

Mz = D(kz + vice) + r 6 + 7 t 6 

in which 


0' 

K‘Z = — 
Oo 


r’P 

K0 = 

r ao 

Ts = J _ y2 f h (e«p + " «p)PC 


(8c) 


(8d) 


The rotationally symmetric external junction loads H* and 
Mz* may arise through the action of an elastic constraint 
such as a reinforcing ring, or they may be independent of the 
shell deformations. Therefore, let 


H e — $n(f€d, c ) + Swfic + H 

Mf = S21 (f C0,c) + S22/3 c + 



and 

5 4 ArH) a + 5 4l 6 (rff)* = (5 4 , £ - , t )(rH) e + 

(5fi, e “ 55,o — 65 t b — fS 2 l5 2 ,e — S22 ) fie + 

(5 6 ,c — 5 6 ,o — 5 6> * — fs 2 i5i, c — Mt) (10d 

in which 

* fwu - r(r'/r)\ * ( v ) 2 1 \ 

*•' = V «c~h Ki ~ \*d), 

•“ ■ (s), 


= ( X1 
' \aoC 


+ X3 + T 3 + 7T3 


5e,j 




Z2 


+ Ts + 7T 5 


* /u) 22 - v(r'/r)\ 

= v «./d A 

), -I! 


(10e) 


Substituting (r//) a and (r//)^ from Eqs. (10a) and (10b) 
into Eqs. (10c) and (10d) yields 


The coefficients Sn to S22 are spring constants, which can be 
determined from the geometry and material properties of the 
elastic constraint. They define the linear dependence of H e 
and M/ on the deformation of shell c; Mz may result from 
nonconcurrence of the three reference meridians. 

When only two shells are involved, Eqs. (7) may be easily 
specialized by deleting the terms with the subscript b. The 
numerical solution of these modified equations will then be 
similar to that of the following section. 


2 . Numerical solution 

Before substituting Eqs. (8) into Eqs. (7), an expression 
relating y/ to y,- is derived. In central difference notation, 
using (4b), 

, = y i+ i — y,_i = (P,_i - Pr'hi + (Pr l <li - 
Vi 2(A© 2(A|) 


or 


y/ = w a , + 


Using Eqs. (4c), 

(Pi - 1 - Pr 1 ) 


w, = 


2(A© 


(J + A,--P,)P,_, - y (9a) 

2(A© 

r = (Pr 1 ?) - g.-i) = _ 

2(A© 

(/ + Aj-'F^qt-i - Aj-'g, 

2(A|) 


and 


7i (rH), + Jzfie = 7s 


7 i(rH) e + 7tfi, = 76 


(Ha) 


(lib) 


Equations (11) will be satisfied independently of (rH), and 
if 

7i = 72 = Ts = 74 = 76 = 7e = 0 


or 

5i, c = 
5*2, 6 = 
5 3 ,c = 


«1.A 


5i,a + 5 i, 6 + f 2 Sn5i. a 5i,b 
5i,a52.6 ~h 5 2 , a 5i,b ~h fsndi,gdi,b 

5l,o + §i,b + f 2 Sinda,§ib 

dl,a&Z,b + 53,q5l,b + fHdigdi'b 

5l,q + Sl,» 4“ f 2 Sll5l,q5l.fc 


+ *«) 

6 s ,e = 5 2 „(|-; + I; + fs 21 ) - 

( 5 2 , q5 4 , q . 5 2tb 5 4l 6 * js \ 

XT + TT “ 56 '“ ■ 

5m = + ^ + 

( 5 3 ,q5 4l q 53,65 4 ,b > 5 |7 \ 

ITT + ■ *') 


( 12 ) 


Expanding the first of Eqs. (9a) gives 

(rH)/ = wn,j(rH)j + Wajfii + x Xti 
/V = Wii, i(rH)i + W 22 ./& + x 2 ,i 


(9b) 


After Eqs. (7a) and (9) are substituted into Eqs. (8), the 


Using Eqs. (12), Eqs. (lOe) may be solved for the com- 
ponents of the W c and x e matrices. Finally, Eqs. (9a) may 
be inverted to give 

Pc-. = {(A + P)-‘[2(A£)AW + B]}.} 

, \ (13) 

2.-1 = — {(A + F)-‘[2(A*)Ax - 2]}. j 



NOVEMBER 1963 


ELASTIC-PLASTIC THIN SHELLS OP REVOLUTION 


2587 


Table 1 Summary of numerical analysis procedure 





Meridional 


Step 

Quantities to be calculated 

stations 

Equation 

1 

J,L,R 



(3) 

2 

r, e, A, T, <!>,¥, Xi to X 4 ,n to r 4 ; xi and xj assumed 
(usually based on and equal to zero) 

[2 to a 
< a -f 2 to6 

(D,(8) 




( c to d — 1 


3 

a,b,f 


Same as above 

(2) 

4 

9 


Same as above 

(2) 

5 

A,,A 4 ,tj to r 5 ;xj to x 

1 

l,a,a *f 1 fiyCyd 

(8) 

6 

kjTnJ 



(3) 

7 

P,q 


fl toa - 1 
\a + 1 toft — 1 

(5) 

8 

WyX 


a and 6 

(9) 

9 

8l,a to 5«,a,5 l,h to 



(10) 

10 

$i,e to 8%, c 

junction calculations 


(12) 

11 

W,x 

P,q \ 


c 

(10) 

12 


c — 1 

(13) 

13 

P,q 


c to d — 1 

(5) 

14 

(rH),, S 


d to c 

(6), (4) 

15 

(rH ),0 — junction calculations 

b and a 

(7), (10) 

16 

(rH),p 


b — 1 to a -h 1 
a - 1 to 1 

(4) 

17 

X*,X 4 


All 

(8) 

18 



All 

(8) 

19 



All 

(14) 

20 



All 

(15) 

21 

XI, X* 


f2 to a 

(1) 


Repeat steps 4 to 19 until convergence of c «p 

(a + 2 to b 



and c£p 


[c to d — 1 


22 



All 

(D,(16) 


Although no station j = c — 1 actually exists in shell c, 
it is convenient to assume its fictitious presence and then to 
use Eqs. (5) to calculate the elements of P c , and the re- 

maining coefficient matrices in the structure. The terms 
(rH) and /3 are then calculated with Eqs. (6) and (4a) for 
shell c, Eqs. (7a, 10a, and 4a) for shell a, and Eqs. (7a, 10b, 
and 4a) for shell 6. Stresses and strains throughout the 
structure may then be calculated from these basic shell 
variables as shown in the following section. Equations (9) 
and (13) are expanded in Ref. 6. 


C. Elastic-Plastic Strain Equations 

Derivation of the following elastic-plastic strain equations 
is given in detail in Ref. 3, in accordance with the von Mises- 
Hencky yield criterion and the deformation theory of plas- 
ticity. After reference surface strains and curvature changes 
have been computed throughout the structure by Eqs. (8), 
strains through the shell walls are calculated by the equations 

ee = €0 + fa € f = ejt + f/q (14a) 

Assuming = 0 and all volume changes are elastic, 


€{• = - i- v ^ + j ZT~ V aT ~ \ ( c fp + € *v) 


(14b) 

By definition, 

e.. = [fe - eo) 2 + (*» - er) 2 + («r - «{) 2 r (14c) 


_ 0_u _ 

- e. E - 


0 


o'. < <r yield 


2(1 + v) <T. ^ . u 

3 <r t > <r yield 


(15a) 


in which the subscript u indicates the results of a uniaxial 
tensile test. A plot of e tp vs may be easily constructed 
using Eqs. (15a). Finally, from the deformation theory of 
plasticity, 


&p = ~ (2ft - et - ft) 


= 3^ (2e* “**-<*) 


(15b) 


Plastic strains computed by Eqs. (15) are used to improve 
the level of approximation in the t terms of Eqs. (lc) and 
(8), establishing the iterative nature of this elastic-plastic 
solution. After convergence of these plastic strains, stresses 
are calculated throughout the structure by Eqs. (Id) and 

<r. = (<r«* — <r$<rs + <r^) Vi (16) 


Results 

A. Summary of Numerical Analysis Procedure 

The method of analysis developed in this report is sum- 
marized in Table 1. This table serves as a convenient flow 
chart for computer programming of the method. The calcu- 
lation procedure is as follows : 

1) Express boundary conditions in finite difference nota- 
tion and determine elements of matrices /, L, and R from 
Eqs. (3). 

2) Calculate geometry terms (r,0,A,T,$, and >P), pressure 
load terms (\ t to X4), and thermal load terms (ri to r 4 ) for 
all interior shell stations, using Eqs. (1) and (8). Assume 
convenient values (usually zero) for plasticity terms and 

3) Calculate elements of matrices A , B, and F for all in- 
terior stations using Eqs. (2). Steps 1 to 3 do not have to 
be repeated during iteration that follows. 

4) Compute elements of g matrix at all interior stations by 
Eqs. (2). 

5) Calculate pressure, thermal, and plasticity terms re- 
quired for strain and moment calculations at boundary and 
junction stations, using Eqs. (8). 

6) Determine elements of remaining boundary condition 
matrices &, m, and t by Eqs. (3). 

7) Calculate elements of coefficient matrices P and q for 
shells a and 6, using Eqs. (5) ; P matrices for these two shells 
do not change during iteration. 

8) To begin junction calculations, determine elements of 
W and x matrices at the junction stations in shells a and b, 
from Eqs. (9). 

9) Calculate d terms for shells a and b using Eqs. (10). 
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a) Geometry and loads. 


2 H 


ELASTIC 

SOLUTION 

— K {»,). K{«„) 

O REFERENCE 7 

ELASTIC-PLASTIC 
SOLUTIONS 
-- K ( ff * fi STRAIN 
K HARDENING 

K (a.)'l PERFECT 
K («.i)J PLASTICITY 



J I I L 



J I L 


J L... 1 .. L_l 


-5 -4 -3 -2 -I 0 -5 -4 -3 -2 -I 0 I 2 3 4 5 

S, DISTANCE FROM JUNCTION. IN. 

b) Effective stress and strain concentration factors. 
Fig. 3 Sample problem. 


10) Calculate 8 terms for shell c using Eqs. (12). 

11) Calculate elements of W and x matrices at junction 
station in shell c by inverting Eqs. (10). 

12) Complete junction calculations by computing elements 
of P and q matrices at fictitious station,; = c — 1. 

13) As a continuation of step 7, calculate elements of P 
and q matrices in shell c. 

14) To begin the calculation of ( rH ) and 0 terms, compute 
their values in shell c, using Eqs. (6) and (4) . 

15) Calculate (rH) and 0 at the junction stations of shells 
a and b by Eqs. (7) and (10). 

16) Continue step 14 in shells a and b. 

17) Calculate plasticity terms 7Ts and n for all shell sta- 
tions using Eqs. (8). 

18) Calculate reference surface strains and curvature 
changes using Eqs. (8). 

19) Calculate total strains throughout the structure by 
Eqs. (14). 

20) Calculate plastic components of strain by Eqs. (15) 
and an equivalent plastic strain vs equivalent total strain 
diagram. 

21) Compute plasticity terms 7Ti and 7 t 2 for all interior 
stations, using Eqs. (1) and the results of step 20. Repeat 
steps 4 to 19 until the plastic components of strain have 
converged sufficiently or for a fixed number of iterations. 

22) Compute stresses throughout the structure, using Eqs. 
(1) and (16). 

For an elastic shell structure, equate all plasticity terms 
to zero and proceed directly from step 19 to step 22. 


B. Sample Problem 

The geometry and loading of the sample problem structure 
are illustrated in Fig. 3. The structure consists of two 


cylinders and a portion of a sphere joined at a common sec- 
tion. Each shell has constant wall thickness and internal 
pressure. The following dimensions, loads, and material 
properties are selected so that yielding is incipient in the 
membrane regions of the structure: f = 50.0 in., h a = 
0.0866 in., fa = 0.0500 in., fa =* 0.1732 in., E = 10 7 psi, 
Pm = 100 psi, Pi = 100 psi, Pc = 200 psi, v - 0.333, and 
oVi.id = 50,000 psi. 

Postyield material behavior is assumed to follow the rela- 
tionship 

f = Cl + Ctfet + Cz(e*t ) 2 €«* > Ci 

e.p { 

1 = 0 < C 4 

The sample problem is solved for both a strain-hardening and 
a perfectly plastic material. For these two cases, the mate- 
rial coefficients become 1) for strain-hardening, Ci = 1.287 X 
10-3, C 2 = -0.8655, Ci ~ 129.5, and C 4 = 4.445 X 10" 3 ; 
and 2) for perfect plasticity, Ci = —4.445 X 10“ 3 , Ct = 1.00, 
Ci = 0, and Ci = 4.445 X 10" 3 . 

The shells are subdivided into meridional and transverse 
stations with (As) = 0.25 in., a = 99, b = 200, c = 202, 
d — 300, and i‘m*x = 7. 

By referring to Fig. 2b and the junction detail of Fig. 3a, 
the external junction moment becomes 

M = 0.1116 V b - 0.1299 V* - -45.75 lb-in./in. 

Finally, if membrane deformation is assumed at all three 
shell boundaries, the boundary condition matrices are 

J = L = R = k = t = 0 


and 


r ( rN( cos<p)„+r 

"52,592 lb" 

L o J 

L 0 . 


The solution to this problem is presented in Fig. 3b in the 
form of effective stress and strain concentration factors. 
These factors are 

K(tr } = = / \ 

(<T i) membrane \ 50,000 psi / max 

jr / s _ (€«t)max _ / ^ \ 

^ ” (^e |) membrane \0.004445 in./in./ max 

The variation of these concentration factors with meridional 
distance s, measured from the junction, is given for the three 
cases considered. Referring to Fig. 3a, negative values of s 
represent distances measured downward from the junction, 
in shells a and 6. 

In Fig. 3b, the solid lines represent the solution for a per- 
fectly elastic material. Stress and strain concentration 
factors are equal with a maximum value of 1.31 occurring at 
the junction in shell b. Circles indicate the results of a 
standard elastic discontinuity solution that uses edge in- 
fluence coefficients. The dashed lines present results for 
the strain-hardening material. Decreases in stress concen- 
tration are accompanied by corresponding increases in strain 
concentration, which reaches a maximum value of 1.44. 
Finally, the dashed-dotted lines show the elastic-perfectly 
plastic solution. This includes the upper limit of strain 
concentration and lower limit of stress concentration for 
various types of postyield behavior. Maximum effective 
strain concentration is 1.89, again occurring at the junction 
in the hemispherical bulkhead. 

The elastic solution is clearly verified by the standard 
discontinuity analysis presented in Ref. 7. The two elastic- 
plastic solutions may be evaluated by determining the ex- 
tent to which they satisfy the junction continuity conditions 
given in Eqs. (7). The quantities that are compared in 
these equations are presented in Table 2. 
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Table 2 Junction deformations and forces 


Quantity 

Strain-hardening solution 

Perfectly plastic solution 


Station 



Station 


a 

b 

c 

a 

b 

c 

0, rad 

0.01439 

0.01439 

0.01439 

0.01671 

0.01671 

0.01671 

eg, in./in. 

0.004522 

0.004573 

0.004540 

0.004795 

0.004819 

0.004814 

(rH\ in. -lb /in. 

-724.42 

536.27 

-188.15 

-404.45 

77.91 

-326.54 

Af, lb-in. /in. 

16.97 

-5.09 

-34.30 

12.22 

-0.93 

-31.75 


The meridional couples were determined from the 
meridional stresses by numerical integration. With 
these quantities, it can be shown that Eqs. (7) are satisfied 
very closely for both cases. 

For further verification, the present method was applied 
to a continuous shell problem for which an elastic-plastic 
solution had already been obtained. This problem, which 
was presented in detail in Ref. 3, consisted of the deformation 
of a uniform cylindrical shell by an axial thermal gradient. 
Excellent agreement between the results of the two methods 
was found, which shows that the numerical analysis pre- 
sented in this report is equally applicable to continuous and 
discontinuous shell regions. 

Conclusions 

A numerical method has been presented for solving prob- 
lems involving the axisymmetric elastic-plastic deformation 
of general shells of revolution which may contain discon- 
tinuities in geometry, material properties, or loads. The 
method was used to determine effective stress and strain 
concentration factors at the common junction of three dis- 
similar shells in the elastic-plastic regime. The following 
conclusions were found : 

1) This method produced an elastic solution that agreed 
very well with that obtained by a standard discontinuity 
analysis. 

2) By using successive approximations, the method pro- 
duced a convergent elastic-plastic solution. 


3) In this elastic-plastic solution, conditions of equilibrium 
and compatibility at the discontinuity were satisfied. 

4) When applied to the analysis of a continuous shell re- 
gion, close agreement with known elastic-plastic results was 
found. 

5) Elastic-plastic shells with discontinuities may be 
analyzed with almost the same directness and ease as com- 
pletely elastic continuous shells. 
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